Persistent cyclonic structures in self-similar turbulent flows 
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Invariance properties of a physical system govern its behavior: energy conservation in turbulence 
drives a wide distribution of energy among modes, as observed in geophysics, astrophysics and 
engineering. In hydrodynamic turbulence, the role of helicity (which measures departures from 
mirror symmetry) remains unclear since it does not alter this distribution. However, the interplay 
of rotation and helicity leads to significant differences. Using numerical simulations we show the 
occurence of long-lived laminar cyclonic vortices together with turbulent vortices, reminiscent of 
recent tornado observations. Furthermore, the small scales are completely self-similar with no 
deviations from Gaussianity. This result points to the discovery of a small parameter in rotating 
helical turbulence. 
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Self-similar behavior in turbulence is viewed as a large- 
scale phenomenon, recently linked to conformal invari- 
ance (e.g., local scale invariance through transformations 
that preserve angles but not distances) when examin- 
ing the scaling properties of vortex lines [l| for two- 
dimensional turbulence Q , and surface quasi-geostrophic 
turbulence Q for a rotating stably stratified layer 
The link found to percolation theory allows for the an- 
alytical determination of scaling exponents such as the 
fractal dimension of vortex clusters, although the con- 
sequences for statistical measures of turbulence, e.g., 
through scaling laws for correlation functions, has not 
been clarified yet. However, the three-dimensional case 
is known to be much more complex. In three dimensional 
turbulence, the flow is not scale invariant and the knowl- 
edge of one exponent does not allow the prediction of the 
exponents for all orders. The search for self-similar quan- 
tities in three-dimensional turbulence is a long-standing 
problem; it would relate its study with critical phenom- 
ena and the out-of-equilibrium statistics of systems with 
a large number of modes, and it would allow the use of 
tools (such as the renormalization group from quan- 
tum field theory and statistical mechanics to further our 
understanding of such flows. 

One of the most important and useful principles of 
physics is that of conservation laws linked, through the 
theorem of Emma Noether, to invariance properties of 
the underlying equations. Conservation of energy is in- 
voked when explaining the observation of a wide range of 
excited scales in a turbulent flow: the nonlinear coupling 
due to advection leads to the feeding of modes at all the 
scales available to the system [|| ; this constant flux direct 
cascade of energy to small scales is arrested by dissipa- 
tive processes. The conservation of helicity H = (u • lj), 
the correlation between the velocity u and the vorticity 
u> = V x u, was discovered much later Q, although it 
appears rather clearly when formulating the Euler equa- 
tions governing the dynamics of the ideal incompressible 
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FIG. 1: Temporal evolution of the maximum vorticity in the 
flow, with uirms ~ 20. The insert gives a short-time blow-up 
in which the inertial wave period 1/n with fi = 9 is slightly 
visible; note the stong local bursts in Lj m ax- 



fluid in terms of the so-called Lamb vector u x u>. Helic- 
ity is not positive definite, unlike energy; it is a topolog- 
ical invariant, representing the degree of knottedness of 
vortex lines [7(, and it is a pseudo-scalar (the symmetry 
group related to its conservation is discussed in Q ) . He- 
licity in turbulent flows can lead to drag reduction, and 
to better mixing of chemical components 0]. 

Turbulent flows are known to develop helical structures 
which arc persistent since their associated non-linear ad- 
vection is weak ; thus, the evolution of these struc- 
tures takes place on the slow dissipative time scale. In the 
atmosphere, the persistence of helical flows was invoked 
to explain structures encountered in supercell storms that 
can give rise to tornadoes 



llj . Even though small-scale 



structures are found to be strongly helical, the helicity 
of a flow does not seem to alter its dynamics: indeed, 
numerical evidence stemming from direct numerical sim- 
ulations of incompressible, isotropic, and homogeneous 
turbulence indicate (both for weak or strong global helic- 
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FIG. 2: Characteristic (integral) length scales for the energy 
(solid line, L) and the helicity (dash line, L H ). The significant 
growth of L confirms the occurrence of a build-up of energy 
at large scale, in what is called an inverse cascade. 

ity) that the distribution of energy among scales follows 
a power-law Q which, expressed in terms of correlation 
functions or of structure functions of second order, reads 
S 2 {1) ~ £ 2 / 3 , with S p (£) = ([u L {r + £) - u L (r)]P) the 
pth-order longitudinal structure function on a distance £, 
ul being the projection of the velocity field along the vec- 
tor £. Similarly, it was shown using the renormalization 
group that the helical contribution to eddy viscosity 
is sub-dominant. 

The question we now address is the interplay of helic- 
ity with rotation, included in the incompressible Navier- 
Stokes equations through the Coriolis force: 

-^+wxu + 2flxu=-VP + yV 2 u + F, (1) 

at 

with V • u = 0; we choose the rotation axis to be in the 
z direction: Jl = Qz, with f2 the rotation frequency; V is 
the total pressure divided by the (constant) mass density 
and modified by the centrifugal term, v is the kinematic 
viscosity and F is an external mechanical force that drives 
the turbulence (mimicking for example a convectivc input 
of energy) ; it is given by the ABC flow [l3| (a Beltrami 
flow with u) proportional to v) with hp = 2-njkp the 
characteristic scale of the forcing, of amplitude F$: 

F = F {[Bcos(k F y) + Cs'm(k F z)]x+ 
+ [C cos(kpz) + A sm(k F x)] y + 
+ [Acos(k F x) + Bsm(k F y)]z} . (2) 

The above equations are integrated with a pseudo- 
spectral code and periodic boundary conditions using a 
second-order Runge-Kutta temporal scheme. Computa- 
tions were performed on a grid of 1536 3 points. The 
flow is first led to establish a statistically steady state 
with = 0.06, i.e. in the near absence of rotation, a 
phase taking roughly ten turn-over times; with k F = 7, 
the rms velocity is U « 1 and the Reynolds number 




FIG. 3: (Color online) Top left: visualization of helicity 
density at t ~ 21; light (green) is negative and dark (blue- 
magenta) is positive. Right: vorticity intensity in the same 
region, with fluid trajectories in red. Note the co-location of 
laminar structures, with smooth paths, and a tangle of vortex 
filaments with more complex paths. Bottom left: A zoom on 
the vorticity magnitude, illustrating the co-existence of both 
structures at very disparate scales. Small-scale vortex fila- 
ments to the left entangle to form filaments at larger scales. 

R v = ULpjv w 5600 with the choice oiv= 1.6 x 10~ 4 . 
Then, at a time labeled t = in the following, the 
rotation is set to fl = 9, corresponding to a Rossby 
number Ro = U/(2ilL F ) w 0.06. The time step is 
At = 2.5 x 10~ 4 . The computation is then performed 
for 30 turn-over times t^l = Lp/U, corresponding to 
180 in units of a large-scale inertial wave period. 

An example of the temporal evolution of the flow is 
given in Fig. [TJ showing the maximum of vorticity, the 
rms value of vorticity being uj rms « 20; a blow-up given 
in the inset shows a slight signature of the period of iner- 
tial waves, with of the order of nine crests per unit time; 
also note that the vorticity in the flow is intermittent, 
both in its time evolution as in its distribution in space. 
Fig. [5] gives as a function of time the integral scale for 
the velocity (solid line) and for the helicity (dash line): 
the monotonous increase in the characteristic scale of the 
energy is a manifestation of a self-organization process of 
the flow at large scales, into structures that contain most 
of the energy. This constant flux transfer of energy to- 
wards large scales is known as an inverse cascade. On 
the other hand, the characteristic scale of the helicity re- 
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mains approximately constant, a manifestation of helicity 
being transferee! towards small scales. 

Two results are striking; they concern the structures 
that develop in the flow and the ensuing statistics. The 
spatial helicity distribution and vorticity magnitude in 
a box of size 0.25 x 0.5 x 1 the total fluid volume, and 
with the vertical direction upward, is given in Fig. [3] we 
observe the juxtaposition in space, for the same physical 
variables, of columnar laminar structures, and of strong 
vortex filaments in a tangled network. Within the column 
there is a strong alignment of v and u) with a resulting 
strong local helicity density. As a result, particle trajec- 
tories (in red) are helical and upward: these are cyclonic 
events that are persistent in time because of the lack of 
nonlincaritics in their vicinity. At any given time, two or 
three such columns can be easily identified in the entire 
domain. The structures are long-lived, and were tracked 
for over ten turnover times without being substantially 
deformed or destroyed by the surrounding turbulent flow. 
On the other hand, the tangle of vortex filaments right 
next to these laminar structures has complex particle tra- 
jectories, lives of the order of one eddy turn-over time and 
corresponds to a classical turbulent fluid albeit with some 
perceptible vertical organization (which is also found in 
the case of non-helical rotating fluids 

Overall, the vorticity is strongest in the tangle of vor- 
tex filaments, whereas the vertical velocity is coherent 
and strongest in the laminar columns (not shown). The 
core of the column is surrounded by a calm region with 
weak vorticity (note the emptiness of vortical structures 
in the surroundings of the column) which acts as a transi- 
tion region between the laminar and the turbulent flow. 
The origin of these stable structures can be identified 
by integrating backward in time the particle trajecto- 
ries: they correspond to regions of large helicity where 
the columns form. Such a laminar organization of the 
velocity field is not observed in helical isotropic and ho- 
mogeneous turbulence, nor is it observed in rotating flows 
without helicity. The interplay between rotation (which 
breaks the mirror-symmetry in the evolution equations) 
and helicity (which quantifies departures from mirror- 
symmetry of the flow) is the driving agent for the forma- 
tion of strong localized and persistent columnar struc- 
tures, even though helicity is mostly transfered towards 
smaller scales [15| and is itself strongly intermittent (see 
below, Fig. 2]). The tangle of vortex filaments surround- 
ing a laminar structure, together with more complex, 
larger and spiraling features, is reminiscent of observa- 
tions of multiple core vortex tornadoes [3| ■ 

The coexistence of the laminar columns with the tur- 
bulent flow has a strong impact on the statistics of the 
velocity field. This can be measured by studying the 
histograms of velocity and helicity differences over a dis- 
tance £ chosen to be in the small scale inertial range, as 
well as by the computation of the exponents £ p of the 
velocity field structure functions S p (i) ~ £*>* previously 
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FIG. 4: (a) Probability density function of increments for the 
velocity (solid line) and helicity (long dash), normalized by 
their variances, averaged over ten turnover times, and taken 
on a distance i = 0.2 (roughly 22% of the forcing scale Lf, 
i.e., in the inertial range of the direct cascade). The velocity is 
close to a Gaussian (dots) , whereas the helicity displays strong 
wings, characteristic of strong localized events, (b) Scaling ex- 
ponents of the velocity structure functions, with error bars. 
Diamonds are for the flow with Q — 0.06, displaying a clas- 
sical three-dimensional intermittent fluid scaling. A straight 
line with £2 = 2/3 (solid) is shown as a reference; the devi- 
ation from the straight line indicates a departure from self- 
similarity. Crosses are for the flow with SI = 9, averaged over 
twelve directions in the plane perpendicular to the rotation 
axis, and over ten snapshots from t = 20 to t = 30. The dot- 
ted straight line indicates self-similar scaling with £2 = 1.42, 
close to the steepest scaling that can be obtained in such a 
flow, as found in [l5| using a simple argument based on a 
direct helicity cascade mediated by inertial waves. 

defined. The increments were taken in twelve directions 
perpendicular to the axis of rotation, and averaged over 
volume and over ten snapshots of the fields, spanning ten 
turnover times. This gives a total of over 4.3 x 10 11 data 
points for each value of the scalar increment I. 

The resulting histograms (Fig. 3]) are strikingly dif- 
ferent: the velocity increments are close to Gaussian, 
whereas for the helicity increments strong wings obtain, 
characteristic of the presence of intense events and inter- 
mittency. Scaling exponents of the velocity up to p = 7 
are shown in Fig. 0Jb). 111 the run with = 0.06, the 
exponents are those observed in turbulence without rota- 
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tion, with for example (3 = 1 , a trace of the conservation 
of energy in the absence of dissipation in those variables. 
The deviation of the exponents from a straight line indi- 
cates the break down of scale invariance in the flow, as 
is well known for three dimensional turbulence: strong 
events are more probable at small scales. In the run 
with £1 = 9, as times grows the system evolves toward a 
different regime: that of self-similarity where the knowl- 
edge of one exponent (say, that of second-order correla- 
tion functions) gives scaling for all higher-order structure 
functions. Henceforth, the scaling of strong gradients is 
completely determined if one of the exponents is known; 
the dotted line (( p ~ 0.71p) represents the data best, in 
agreement with predictions made in [lH for £2- 

This is a remarkable property, to our knowledge never 
before observed in small-scale three dimensional turbu- 
lence with a direct cascade of energy. While the energy 
has organized in fully helical updrafts at large scales, the 
helicity itself concentrates in the cores of these struc- 
tures as well as in the small-scale tangle of vortex fila- 
ments. The identification of self-similarity in turbulence 
with both helicity and rotation (conditions that are rel- 
evant in many atmospheric flows) would allow us to re- 
late the dynamics of such three-dimensional flows to the 
advances made in two-dimensional turbulence and crit- 
ical phenomena in general However, in order to be 
able to use renormalization group techniques 0, and 
other tools from statistical mechanics, a small param- 
eter needs to be identified, besides the Rossby number 
that governs the energetic exchanges between turbulent 
eddies and waves when dealing, e.g., with the weak tur- 
bulence regime 17|. Indeed, the smallncss of the ratio 



of the inertial wave period to the eddy turn-over time 
has already been used to derive intcgro-differcntial equa- 
tions in terms of energy and helicity spectra fl~8| in the 
context of weak turbulence. However, these solutions are 
not observed in our study or in many atmospheric flows 
for at least three reasons: first, the Rossby number in the 
atmosphere and in our simulations is moderate and far 
from the limit considered in the theory; second, the the- 
ory is non-uniform in scale and the weak turbulence limit 
breaks down; third, in the case of rotation, the inverse 
cascade of energy is not present at lowest order in the 
theory and thus the solution selected by this approach is 
one of an energy cascade to small scales, whereas our nu- 
merical data indicates that this cascade is sub-dominant 



to the helicity cascade [15 1. 

This sub-dominant direct energy cascade to the small 
scales provides in fact the needed small parameter for tur- 
bulence with rotation in the presence of helicity, in the 
form of the (adimcnsionalized) ratio \ = e l Lf^, where 
e and e are the energy and helicity fluxes, constant by 
definition in the inertial range. The energy flux to small 
scales is all the more negligible as more energy is trans- 
ferred to large scales in an inverse cascade. The con- 
jecture (now backed by numerical data) that x < 1 as 



the Rossby number gets small (but not necessarily tend- 
ing itself to zero, see 15|, the transition taking place 
for Ro ~ 0.1), led us in fact to the hypothesis that for 
moderate to strong rotation, the small-scale cascade is 
dominated by helicity and mediated by waves, leading to 
a prediction for the scaling £2 = 1.5 for rotating flows 
with maximum helicity 15| . The lack of intermittency 
in the direct energy cascade is thus a novel feature pro- 
viding for the first time a small parameter for small-scale 
three-dimensional turbulent flows, and giving a new way 
to look into turbulence. 
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